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ABSTRACT 

We explore the formation of a-enhanced and metal-rich stellar populations in the nuclei 
of luminous ellipticals under the assumption of two extreme galaxy formation scenarios 
based on hierarchical clustering, namely a fast clumpy collapse and the merger of two 
spirals. We investigate the parameter space of star formation time-scale, IMF slope, 
and stellar yields. In particular, the latter add a huge uncertainty in constraining 
time-s cales and IMF slopes. We find that - for Thielemann, Nomoto & Hashimoto 
( 1996 ) nucleosynthesis - in a fast clumpy collapse scenario an [a/Fe] overabundance 
of ~ 0.2 dex in the high metallicity stars can be achieved with a Salpeter IMF and 
star formation time-scales of the order 10 9 yr. The scenario of two merging spirals 
which are similar to our Galaxy, instead, fails to reproduce a-enhanced abundance 
ratios in the metal-rich stars, unless the IMF is flattened during the burst ignited 
by the merger. This result is independent of the burst time-scale. We suggest that 
abundance gradients give hints to distinguish between the two extreme formation 
scenarios considered in this paper. 

Key words: galaxies: elliptical, formation scenario - a-enhancement - stellar popu- 
lations - initial mass function 



1 INTRODUCTION 



An enhancement of a-elements^| in the metal-poor stars of 
our Galaxy, with respect to the s olar p roportions, wa s firs t 
realized by Alter fe G reenstein ( |l96o| , Wallerstein ( |1962| ) 



and Conti et al. ( 1967 ). Subsequent studies (see Truran & 
Burkert 1993; Gehren 1995, and references therein) have 
confirmed that the abundance ratios in metal-poor stars ex- 
hibit a/Fe larger than the solar ratios. Since a-elements are 
mainly produced in short-lived, massive stars, while a sub- 
stantial contribution to Fe comes from Type la supernovae 
(SN la) with longer lived progenitors, this pattern is gen- 
erally interpreted in terms of formation time-scales. In the 
early stages of the chemical evolution (i.e. at low metallic- 
ities) the stellar abundance ratios reflect the SN II prod- 
ucts. As the metallicity in the interstellar medium (ISM) 
builds up towards the solar value, the [Mg/Fe] ratio is sig- 
nificantly driven down by the products of SN la explosions 
(cf. McWilliam 1997 and references therein). Abundance ra- 
tios therefore carry information about the main contributors 
to the enrichment and the star formation time-scales (SFT). 
In contrast to the Milky Way, it is generally impossi- 



* The so-called a-elements, i.e. 
by synthesizing 4 He-particles. 



O, Mg, Si, Ca, Ti, are build up 



ble to resolve stars in elliptical galaxies, thus their chemical 
properties have to be inferred from spectral indices. The- 
oretical population synthesis models based on solar abun- 
dance ratios, while adequately describing the (Fe)-Mg2 fea- 
tures observed in the low luminosity ellipticals, fail to repro- 
duce these indices in brighter objects (Peletier 1989; Gorgas, 
Efstathiou & Aragon Salamanca 1990; Faber, Worthey & 
Gonzalez 1992; Worthey, Faber & Gonzalez 1992; Davies, 
Sadler & Peletier 1993). Specifically, at a given (Fe) index 
the Mg2 index is stronger than predicted by such models. 
This result is usually interpreted in terms of an [Mg/Fe] 
overabundance in the stellar populations inhabiting the 
most massive ellipticals. The following explanations for this 
Q-enhancement have been proposed (e.g. Faber et al. 1992): 

(i) Short time-scales characterizing the process of forma- 
tion of the stellar populations (SP) dominating the observed 
line indices. 

(ii) Selective mass loss, such that SN la products are lost 
more efficiently. 

(iii) An Initial Mass Function (IMF) biased towards mas- 
sive stars. 

(iv) Lower rate of SN la with respect to what is appro- 
priate for the Milky Way. 

In cases (i) and (ii) the intracluster medium (ICM) should 
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be overabundant in iron, if the contribution of giant ellipti- 
cals (h ereaf ter gEs) to the ICM pollution is large (Renzini 
et al. 1993). Recent results from ASCA indicate that the 



element ratios in the ICM are consistent with solar values 
(Mushotzky et al. 1996; Ishimaru & Arimoto 1997; Gibson, 
Loewenstein & Mushotzky 1997; Renzini 1997). This sug- 
gests that either gEs do not play a major role in establishing 
the ICM composition or options (i) and (ii) do not apply. 
Case (iv) is an ad hoc assumption without real justification. 
Observed nova rates and LMXB populations make it un- 
likely that early-type ga laxies form less binary stars than 



later types ( Davies 199£ ) 



by several authors ( 


Matteucci 1994; Tantalo et al. 19961 


Vazdekis et al. 1996; Vazdekis et al. 1997; |Gibson & Mat- 


teucci 4997 


; Chiosi et al. 1997J) in the framework of mono- 



lithic collapse models for the galaxy formation. These mod- 
els describe the elliptical galaxy as a single body, and their 
results are hence appropriate to describe the global prop- 
erties of these objects. These studies concentrate on repro- 
ducing the trend of a-enhancement increasing with galactic 
mass, that can be inferred from the observed spectral in- 
dices. The main results of these computations are: 

(i) The classical wind models yield decreasing [Mg/Fe] 



overabi ndance with increasing galactic mass in contradic- 



tion to the observations. 

(ii) An increasing [Mg/Fe] with mass can be reproduced 
if the SFT becomes shorter with increasing mass of the ellip- 
tical. This is t he cas e for the inverse wind models proposed 
by Matteucci (|l994|). 



(iii) Alternatively, an IMF getting flatter with increasing 
galactic mass reproduces the same pattern. 

The most recent works o n this subject tend to favour the 
flat IMF option for gEs (|Vazdekis et al. 1997j ; |Gibson & 
Matteucci 1997j ; |Chiosi et al. 1997j ). 

In this paper we re-investigate options (i) and (iii). 
Rather than addressing the question of the increasing over- 
abund ance with galactic mass, we seek under which condi- 



tions ai overabundance at high Z - appropriate for nuclei of 
gEs - is produced. In doing so we consider a physically plau- 
sible scenario for the formation of gian t ellipticals based on 



hierarchical cluste ring of CDM theory ([White fc Rees 1978 



Frenk et al. 1985| ; |Efstathiou et al. 1988[ ). 

Indeed there are several indications that merging must 
play an important role in galaxy formation. Approximately 
50 per cent of luminous e l lipticals host kinematically decou- 
pled cores (Bender 199C; Bender 1996). In some cases, the 



SPs in these cor es show markedly differ e nt properties in th e 



spectral indices (Bender & Surma 1992; Davies et al. 1993) 



As shown by Greggio (1997, hereafter G97), the strong Mg2 
and (F e) indices measured in the nuclei of gEs indicate a 



substar tial degree of pre-enrichment in the gas which is con- 



verted into stars in the inner parts of these galaxies. Al- 
though not exclusively, this can be accomplished in a merg- 
ing formation scenario, provided that gas dissipation plays 
an important role. 

Ki nematical properties and line strength features are 



then de termined by the ratio of star formation and merg- 
ing time-scales (Bender & Surma 1992; Bender, Burstein & 
Faber 1992) . In this framework we can envisage two ex- 
treme cases for the formation scenario for ellipticals, mainly 



distinguished by the overall formation time-scale of the ob- 
ject: 

a) Fast Clumpy Collapse. On rather short time-scales (~ 1 
Gyr), massive objects are built up by merging of smaller 
entities. Star formation (SF) occurs within these entities as 
they merge. The different dissipative properties of the newly 
formed stars and the gas produce a chemical separation, 
such that the enriched material flows to the centre (Bender 
& Surma 1992; G97). On a larger scale, however, the whole 
system participates in a general collapse, merging and SF 
occuring within the same (short) time-scale. We notice that, 
from a chemical evolution point of view, this formation mode 
can also be considered as a refinement of the classical Larson 



(1974) monolithic collapse, 
b) Merging Spirals. The merging event occurs when the 
merging entities have already converted most of their gas 
into stars. For example, two spiral s similar the Milky Way 
merge to form an elliptical gal axy ([Farouki fc Shapiro 198S ; 
Ncgropontc fc White 1983] ; perhard 1983| ; |Barnes""l98£ ; 
Hernquist 1993[ ). Most of the stars in the merging entities 
have formed in a continuous and long lasting (~ 10 Gyr) SF 
process, leading to approximately solar abundances in the 
ISM. At me rging, the (enriched) resid u al gas flows down to 
the c entre (Hernquist & Barnes 1991; Barnes & Hernquist 



199f ) where it experiences a violent SF episode, with a short 
time-scale. The global formation process, however, lasts ~ 
10 Gyr or more. 

Semi-analytic models of galaxy formation in the framework 
of the CDM theory for structure formation (Kauffmann, 
White & Guiderdoni 1993; Lacey et al. 1993; Cole et al. 
1994) provide a continuum between these two extremes, 
with a distribution of time-scales for the formation of el- 
lipticals which depends on specific assumptions in the mod- 
eling. There are several arguments constraining short forma- 
tion time-scales and old formation ages of elliptical galaxies, 
hence pointing towards the fast clumpy collapse model. They 
are mainly based on the tightness of Fundam ental Plane re- 
lations, i.e. the relation between SPs and o (Dressier et al 
1987; Djorgovski & Davis 1987) an d the small scatter i n 
M/ L ratios perpendicular t o the FP (Renzini fc Ciotti 1993 ). 
Also the colour evolution (Aragon Salamanca et al. 1993), 



and the evolution of the Mg-cr relation with redshift (Ben- 
der, Ziegler & Bruzual 1996; Ziegler & Bender 1997) put 
the star formation ages of cluster ellipticals to high red- 
shift (z J> 2). Despite this, the chemical outcome of sce- 
nario (b) is of interest, since we know that such mergers 
do occur (IRAS galaxies, e.g. Soifer et al. 1984; Joseph & 
Wright 1985; Sanders et al. 1988; Melnick & Mirabel 1990). 
If ellipticals in the field and in clusters are intrinsically dif- 
ferent (e.g. Kauffmann, Chariot & White 1996), an d those 



in clusters form earlier, and on shorter time-scales (Kauff- 



mann 1996), the two distinct scenarios may be assigned to 
these two kinds of objects. We notice, however, that the a- 
cnhancement does not seem to depend on cluster properties 
(J0rgensen, Franx & Kjaergaard 1995), at least within the 
errors. Also, the Mg-cr relation measured in many cl usters is 
remarkably independent of the clusters properties (Colless 



et al. 1998 ). Finally, the Mg-cr relation for field ellipticals 
does not differ appreciably from that derived for the Coma 
Cluster (Bernardi et al., in preparation). Despite this, the 
observational situation is still far from being settled, and 
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we investigate on the predictions of both models to help in 
constraining the problem of the formation of ellipticals. 

The paper is organized as follows. In Section ^ we pro- 
vide a more detailed prescription of the theoretical model 
used in this paper. The results are presented and discussed 
in Section |^ and Section ^, respectively. The main conclu- 
sions are summarized in Section 



2 THE THEORETICAL MODEL 

The basic constraints on our modeling, aimed to describe 
the SPs in the central parts of bright ellipticals, are: 

(i) Achieving high (super-solar) total metallicities (Z); 

(ii) Producing an Q-enhancement of the order of [a/Fe]~ 
0.2-0.4 dex, as implied by the observations. 

The first requirement comes from the observed nuclear Mg2 
index and points towards either higher yields (i.e. shallow 
IMF) in the central regions or to a scenario of enriched inflow 
(Edmunds 1990, 1992; G97). The latter implies that the SPs 
inhabiting the nuclear regions of ellipticals are characterized 
by a metallicity distribution with a minimum metallicity 
Zm > (i.e. they are composite stellar populations, here- 
inafter CSP). The second requirement, instead, comes from 
both Mg and Fe indices, thus in principle the two constraints 
might refer to different components of the CSP. However, 
since the Mg and Fe indices are measured in the same spec- 
tral range, the component of the CSP dominating the Mg2 
index also dominates the Fe index. Hence, requirements 
(i) and (ii) must be met simultaneously from the SP 
dominating the visual light. 

The first constraint is accomplished with a high degree 
of chemical processing, hence a large fractional amount of 
gas has to be turned into stellar mass. This can be achieved 
either in short time-scales (perhaps like in galactic bulges) 
or long time-scales (as in galactic disks). The second char- 
acteristic points towards relatively short time-scales. Note 
that both constraints are generally more easily met when 
assuming a flatter IMF. 

The fast clumpy collapse and the merging spirals modes 
differ mainly in their time-scales of the pre-enrichment pro- 
cesses. These are mimicked by assuming different initial con- 
ditions for the chemical processing of the gas which forms 
the CSP hosted in the central parts of the galaxy. In the fast 
clumpy collapse the process of formation is so fast that the 
galaxy can be viewed as a single body. Its chemical evolution 
can be described as a closed box with only radial separation 
of different SPs. The initial abundance for the SF process is 
assumed to be primordial. In the merging spirals model, in- 
stead, the gas which is turned into stars at the merging event 
has been enriched in the extended evolution history of the 
parent spirals leading to approximately solar abundances. 
The two different star formation histories are sketched in 
Fig. |l| The duration of the active SF process will be ad- 
dressed to as the star formation time-scale tsp (arrows in 
the figure). In the case of the merging spirals model (bottom 
panel) the overall formation time-scale for the final elliptical 
is much longer than tsp. The phase of chemical enrichment 
in the parent spirals before the merger is indicated by the 
dashed line. The chemical outcome of this phase is reflected 
in the initial abundances assumed for the gas participating 




FAST CLUMPY COLLAPSE 



MERGING SPIRALS 



overall time — scale 



evolutionary time 

Figure 1. A sketch of the star formation histories assumed for 
the fast clumpy collapse (top panel) and the merging spirals (bot- 
tom panel) schemes. The arrows indicate the duration of the SF 
episode in the modeling (rgp). The thick solid lines indicate the 
SF episode yielding the SPs which end up in the central parts of 
the modeled galaxy. The dashed line (bottom panel only) sketches 
the rough star formation history of the parent spirals before the 
merger. 

the merger (see Table §). The thick lines in Fig. | indicate 
the metal-rich SPs that are assumed to form the nucleus of 
the object. 

In the simulations, we consider a normalized mass of 
gas and calculate its chemical evolution until the gas is com- 
pletely used up, yielding the maximum metallicity, which is 
appropriate for the most massive ellipticals. It shoul d be 
noticed that in the mode ls by Arimoto & Yoshii (1987) and 
Matteucci & Tornambe ( 1987 ) at the onset of the galactic 
wind the gas fraction is 0.01 — 0.05 in the most massive 
models, depending on the IMF slope. 

One further distinction between the two extreme cases 
for elliptical galaxy formation concerns the contribution of 
SNe la to the Fe enrichment of the high Z stars. In the 
merger of the two spirals we neglect the contribution from 
those SNe la whose progenitors were born in the parent spi- 
rals (dashed line in Fig. |l]), and which explode after the 
merging. In other words, we assume that the products of 
these SNe la are not efficiently mixed down to the central 
part where the final burst occurs. This corresponds to max- 
imizing the [Mg/Fe] overabundance achieved in the nuclear 
CSP. In a fast clumpy collapse, instead, the merging enti- 
ties are in close interaction, chemical mixing is likely to be 
more efficient. The model corresponds to the formation of 
the whole elliptical in a single event, and it seems quite ar- 
tificial to exclude a fraction of the SN la products from the 
enrichment. 



2.1 Input ingredients for the chemical evolution 

We adopt a constant SFR in our computations, since in both 
models the SF episode generating the high Z stars has to be 
short to accomplish significant [Mg/Fe] overabundances. 

The chemical evolution is calculated by solving the 
usual set of differential equations as described in Thomas, 
Greggio fe Bender (1998, hereafter Paper I). Enrichment due 
to planetary nebulae from low-mass stars (Rcnzini fe Voli 
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1981), SNe II from high- mass stars (Woosley & Weaver 1995, 
hereafter WW95, and Thielemann, Nomoto & Hashimoto 
1996, hereafter TNH96) and SNe la occuring in close binary 
systems (Nomoto, Thielemann & Yokoi 1984) are taken into 
account. The fraction A of close binaries is calibrated on a 
simultaneous fit of the age-metallicity relation in the solar 
neighbourhood and the ratio Nin/Nu ~ 5 typically observed 
in Sbc galaxies (van den Bergh & Tammann 1991), leading 
to A = 0.035 (see Paper I). The IMF is assumed as the 
usual single slope power-law normalized by mass, hence th e 
slope x = 1.35 refers to the Salpeter value ( Salpeter 1955 ). 
The upper and lower mass cut-offs are 40 Mq and 0.1 M®, 
respectively. The effect of adopting a different contribution 
from SN la, different mass cut-offs, and an IMF deviating 
from the single power law (e.g. Scalo 1986) are discussed 
below. 

The time-scale tsf, and the IMF slope x represent the 
main parameters in the calculations. 



2.2 Average abundances 

The output of our models are abundances. For a compar- 
ison to the observational data the luminosity weighted av- 
erage abundances have to be derived. The real abundances 
averaged by mass are larger, since at constant age metal- 
poor stars are brighter (e.g. G97). We shall compute both 
averages to figure out the effect quantitatively. 

The mass averaged abundance of element i in a CSP of 
age t g , including living stars and remnants, is given by the 
following eqn.: 



f£Xi(t)ip(t) m s t SP (t g -t)dt 
r n ^(t) mf sp (t B -t)dt 



(1) 



where to and ti are respectively the epochs of start and 
end of the active SF process. The quantity mf p (( 9 — t) 
is the fraction of mass in stars (living and remnants) that, 
being born in t, are present at the epoch t g (SSP means 
Simple Stellar Population of a fixed single age and single 
metallicity): 



ml SP (t g -t) 

tnta(t g —t) 



(m) dm + 



(m) dm 



(2) 



"Ho (tg—t) 



In Eqn. y mto is the turn-off mass at epoch t g of the SSP 
born at time t, m m i n and m max are respectively the lower and 
upper mass cut-offs of the IMF, and w m is the remnant mass 
of a star of initial mass m. The values of function of 

mass are taken from the nucleosynthesis prescriptions (see 
Paper I). As usual, the IMF 4>(m) is normalized to 1. Thus 
the product ip(t) x mf SP (t g — t) gives the mass contributed 
at the time t g from the stellar generation born at epoch t, 
when the abundance in the ISM was Xi(t). 

In the light averaged abundance, the weights are given 
by the luminosity contributed at epoch t g from the SSP born 
at time t: 



(Xi) v (t g ) = 



f^Xi(t)tp(t) Lf p (t g ~t)dt 
Lf p (t g -t)dt 



(3) 



where we have specified the V-band wavelength range be- 
cause this is most relevant to our application. L^p p (t g — t) is 



Table 1. The initial abundances of the gas assumed to form 
the central stellar populations of the elliptical according to the 
merger of two spirals. The values for 0.5 Zq are adopted from 
Paper I. The quantity r prc denotes the evolutionary time-scale of 
the parent spirals before the merger. 



Z[ n Tprc 

0.5 Z Q 3 Gyr 
1.0 Z Q 10 Gyr 



[Fe/H] in [Mg/Fe] in 
-0.39 0.08 
0.00 0.00 



the luminosity at epoch t g of an SSP of unitary mass with 
an age of (t g — t), and is given by SSP photometric models 
as a function of age and metallic ity. In our application we 
will use the models by Worthey (1994). The spread in age 
and metallicity of the CSP are [to,ti] and [Z m (to), Zm(£i)], 
respectively. 

The averaged abundance ratio is obtained by taking the 
ratio of the averaged abundances. For the sake of clarity, 
however, we use the following notation: 



([Xi/Xj]) EE [(X^/iX,)] 



(4) 



where the square brackets denote the usual normalization 
to the solar abundances, for which we use the meteoritic 



abundances given by Anders & Grevesse ( 198£ ) 



3 RESULTS 

In this section we present the results of our computations 
for the two extreme models of elliptical galaxy formation. In 
particular, we consider the average abundances 15 Gyr after 
the first episode of SF. However, since the passive evolution 
of SPs is most rapid in the first Gyr, the result is insensitive 
t o the details of the star formation history and the exact age 
( |Edmunds 1992[ ). 

Since most of the data in the literature concern mag- 
nesium and iron indices, we will concentrate on the average 
abundances of these two elements, thus on the possibility 
of producing an [Mg/Fe] overabundance. Mg can be consid- 
ered as a representative for a-elements. Indeed, the ratios 
[Mg/Fe] and [O/Fe] do not differ significantly in the TNH96 
models (Paper I). Also, we consider V-luminosity averaged 
quantities, since both the Mg and Fe indices fall in this spec- 
tral range. 

In the following, we first discuss the average abundances 
for the closed box model. This is a good description of the 
global properties of an elliptical galaxy formed in a fast 
clumpy collapse mode. Second we look at the abundances 
of the high metallicity component of the closed box model. 
This is assumed to describe the properties of the nuclear 
regions of an elliptical galaxy formed in a fast clumpy col- 
lapse mode. Third we consider the merging spirals model in 
two possible initial conditions (see Table [j]), corresponding 
to different degrees of enrichment of the ISM at the merging 
event, hence to different pre-enrichment time-scales. 

The complete set of numbers for different parameter 
pairs of tsf, x, and Z m and for both formation scenarios is 
presented in Table B. 
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Table 2. The respective average stellar abundances (V-luminosity weighted) resulting from the fast clumpy collapse 
model (columns 4 — 7) and the merging spirals scheme (columns 8 — 10) for the different parameters IMF slope x, star 
formation time-scale, and minimum metallicity Z m . TNH96 yields are adopted. In column 4, the quantity mcsp denotes 
the mass fraction of the segregated composite stellar population for the fast clumpy collapse. In the merging spiral model, 
these fractions are determined by the fractions of gaseous mass when the merger occurs. The model for the evolution of 
the solar neighbourhood (Paper I) yields 0.5 and 0.2 for Z m = 0.5 Zq and Z m = 1 Zq, respectively. 



FAST CLUMPY COLLAPSE MERGING SPIRALS 



1.35 



1.00 



0.70 



r (Gyr) 


Z m {Zq) 


m CSP 


([Mg/Fe]) v 


([Fc/H]) v 


([Z]>V 


([Mg/Fe]) v 


([Fe/H]) v 




0.3 


0.0 


1.0 


0.27 


-0.52 


-0.27 


_ 


_ 


_ 




0.5 


0.5 


0.26 


-0.20 


0.03 


0.18 


-0.12 


0.05 




1.0 


0.2 


0.25 


-0.03 


0.19 


0.09 


0.14 


0.20 


1.0 


0.0 


1.0 


0.20 


-0.42 


-0.22 


_ 


_ 


_ 




0.5 


0.5 


0.19 


-0.12 


0.06 


0.15 


-0.07 


0.07 




1.0 


0.3 


0.18 


0.03 


0.20 


0.07 


0.18 


0.23 


4.0 


0.0 


1.0 


0.09 


-0.24 


-0.14 


_ 


_ 


_ 




0.5 


0.6 


0.07 


-0.01 


0.09 


0.08 


0.04 


0.12 




1.0 


0.3 


0.06 


0.17 


0.22 


0.04 


0.23 


0.29 


0.3 


0.0 


1.0 


0.34 


-0.05 


0.24 










0.5 


0.8 


0.34 


0.11 


0.38 


0.28 


0.17 


0.39 




1.0 


0.7 


0.33 


0.20 


0.47 


0.21 


0.34 


0.48 


1.0 


0.0 


1.0 


0.28 


0.07 


0.31 










0.5 


0.8 


0.28 


0.20 


0.43 


0.24 


0.25 


0.44 




1.0 


0.7 


0.28 


0.29 


0.50 


0.19 


0.41 


0.53 


4.0 


0.0 


1.0 


0.20 


0.23 


0.39 










0.5 


0.8 


0.20 


0.34 


0.48 


0.18 


0.38 


0.50 




1.0 


0.7 


0.20 


0.42 


0.55 


0.15 


0.55 


0.61 


0.3 


0.0 


1.0 


0.40 


0.31 


0.62 










0.5 


0.9 


0.40 


0.40 


0.70 


0.36 


0.44 


0.70 




1.0 


0.8 


0.40 


0.45 


0.74 


0.32 


0.55 


0.75 


1.0 


0.0 


1.0 


0.35 


0.45 


0.70 










0.5 


0.9 


0.35 


0.52 


0.76 


0.33 


0.55 


0.76 




1.0 


0.9 


0.35 


0.56 


0.79 


0.29 


0.65 


0.81 


4.0 


0.0 


1.0 


0.29 


0.62 


0.79 










0.5 


0.9 


0.29 


0.67 


0.83 


0.27 


0.70 


0.84 




1.0 


0.9 


0.29 


0.71 


0.86 


0.25 


0.84 


0.92 



3.1 Predictions from the closed box model 

As mentioned in the Introduction, an [Mg/Fe] overabun- 
dance can be realized either with short SFTs or with a flat 
IMF. Correspondingly, different IMF slopes constrain dif- 
ferently the tsf parameter. We investigate quantitatively 
this tsf-x degeneracy, and consider two prescriptions for 
the SN II nucleosynthesis (WW95 and TNH96). In case of 
WW95 we adopt the results from model B which assumes an 
enhancement of the explosion energy in high mass stars by a 
factor 1.5 and therefore compares best to the TNH96 mod- 
els. It should be noted that model B produces the highest 
Mg/Fe ratio among the WW95 calculations (models A,B,C). 
For a detailed discussion we refer the reader to Paper I. 



3.1.1 WW95 yields 

Fig. [2| shows the averaged cv-enhancement as a function of 
SFT and IMF slope for the WW95 stellar yields. The pa- 
rameter range for tsf and x are chosen such that the result- 
ing abundance ratios fall in the range indicated by the data 
(0.2 — 0.4 dex). As shown in Paper I, for these stellar yields 
the [Mg/Fe] ratio in the SN II ejecta from one SSP is only 
slightly above solar. It follows that extremely short tsf are 
required to obtain values significantly above solar. In fact, 
Fig. ^ shows that a tsf ~ 10 7 yr is necessary, irrespective of 



the IMF slope. These results are presumably not compatible 
with the time-scale to cool the ejecta and to include them 
in the next generation of stars. 

The minimum in ( [Mg/Fe] )v at r SF = 11-17 Myr is 
a consequence of the metallicity dependence of the WW95 
yields (included in our computations), in which the highest 
[Mg/Fe] ratios are achieved for lowest (10 -4 Zq) and high- 
est (1 Zq) initial metallicity (see tables B1-B12 in Paper I). 
Thus in a regime where SPs of intermediate Z dominate the 
enrichment, the calculated overabundance reaches a mini- 
mum. 

It should be noticed that the model for the SN la rate 
adopted in the computations does not influence this result. 
The Fe contribution from these objects starts to play a role 
only at ages larger than 10 s yr. Lower Fe yields from SN II 
would lead to larger ( [Mg/Fe] )v for a given tsf- However, 
most of the Fe comes from stars of mass ^ 20 Mq for which 
a lower Fe contribution would be unlikely (see Paper I). 



In the model by Matteucci (1994) an a-element over- 
abundance is still produced with SFTs of the order 10 8 yr. 
I t is p ossible that the different re sults come from Matteucci 
( 1994 ) having adopted Woosley ( 1986 ) yields which do not 
account for the fall back in high mass stars, and are there- 
fore characterized by a huge magnesium yield from massive 
stars. We have shown in Paper I (Table 6) that an extrap- 
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Figure 2. The figure shows the average [Mg/Fe] overabundance 
in the stars weighted by V-band luminosity as a function of star 
formation time-scale and IMF slope (by mass). SN II nucleosyn- 
thesis is adopted from WW95. The contour lines mark the respec- 
tive levels of constant ([Mg/Fe])y No segregation of metal-rich 
populations is considered, therefore the results are interpreted 
as global properties of the galaxy. The time-scales required for a 
significant ^-enhancement are unplausiblc short (~ 10 7 yr). 



olation of the WW95 yields to the high mass end does not 
increase the Mg abundance in the ejecta of one SP signifi- 
cantly, again because of the fall back effect. 



3.1.2 TNH96 yields 

The results for TNH96 nucleosynthesis are presented in 
Fig. |^. In this case, tsf as long as a few Gyr still leads 
to an [Mg/Fe] overabundance, provided that the IMF is not 
very steep. It is important to emphasize that the higher 
Mg yields in TNH96 basically stem from a higher contribu- 
tion by stars in the intermediate mass range (18 — 25 Mg) 
which are not severely affected by fall back (Paper I). The 
results are more clearly illustrated in the 2-dimensional pan- 
els of Fig. ^| where we plot contours of constant ([Mg/Fe] )v, 
([Fe/H])v, and {[Z])v in the rsF-x-plane. 

The average Z and [Fe/H] are little sensitive to tsf, 
while they depend on the IMF slope. This is the mere con- 
sequence of having considered a closed box model which uses 
up all of the gas. In these cases the metallicity distributions 
extend up to the maximum possible value, which is given by 
the stellar yields of the SSP, and therefore only depends on 
the IMF slope, for a given stellar nucleosynthesis prescrip- 
tion. The contours are not totally flat because of the relax- 
ation of the instantaneous recycling approximation which 
allows a delay in the enrichment. Since a substantial frac- 
tion of Fe comes from SN la the time dependence is larger for 
([Fe/H])v- The plots also show that the iron abundances are 
lower than the actual total metallicity by roughly 0.1 — 0.2 
dex. Hence at SFTs of the order one Gyr and below, Fe/H 
does not trace total metallicity Z very well. 

From Fig. ^|one can see that ([Fe/H])v ~ is achieved 
with x in the range 1.2 — 0.9. For these values of the IMF 
slope an [Mg/Fe] overabundance of 0.25 is achieved with 
tsf in the range 0.8 to 2.6 Gyr. Restricting to the Salpeter's 
IMF (x = 1.35), one needs tsf ~ 0.1 — 1 Gyr for an a- 



enhancement in the range ( [Mg/Fe] )v ~ 0.2 — 0.3 dex. The 
average iron abundance and metallicity, however, are low 
reaching ([Fe/H]) v ~ -0.4 and ([Z]) v ~ -0.2 for tsf ~ 1 
Gyr. The fact that the Salpeter IMF slope predicts sub-solar 
average metallicities should not be regarded as demanding 
flat IMFs for ellipticals. The observational constraints for 



super-solar metallicities (Davies et al. 1987; Wbrthey et al 



1992; ) concern the central parts of the galaxies, while the 



results discussed here refer to the global average metallicity 
of the galaxy. We will concentrate on this problem in the 
following subsections. 

From now on we restrict our simulations to TNH96 nu- 
cleosynthesis, since, when adopting WW95 yields, the SF 
time-scales required to accomplish a [Mg/Fe] overabundance 
seem unrealistically short. 



3.2 Central abundances for the Clumpy Collapse 
model 

In this model, the stellar populations dominating the central 
light of elliptical galaxies are those in the high metallicity tail 
of the closed box distribution. Since the chemical trajectory 
is characterized by decreasing overabundance as the metal- 
licity increases, the larger the minimum metallicity (Zm) 
adopted for the CSP inhabiting the central parts, the lower 
the average overabundance. This is demonstrated quanti- 
tatively in Fig. ^, where we show the predictions for our 
family of closed box models with Salpeter IMF (left-hand 
panel) and flat IMF with x = 0.7 (right-hand panel), and 
varying tsf. The dotted and dashed lines bracket the pos- 
sible solutions, showing respectively the results for Z m — 
(no segregation of the high Z populations) and Zm = Zm, 
which is the highest Z attained by the closed box model 
(Zm ~ 3 Zq, 17 Zq for x — 1.35 and x — 0.7, respec- 
tively). Adopting the Salpeter slope and tsf ^5 1 Gyr a 
([Mg/Fe]) v ^ 0.2 dex is obtained for every choice of Z m . 
The amount of iron locked in stars, instead, is much more 
sensitive to the segregation of the populations. In the case of 
tsf = 1 Gyr when cutting the low-metallicity tail at 0.5 Zq, 
([Fe/H])v increases from —0.44 to —0.15 dex. For the 
super-solar populations in the centres of ellipticals, 
solar iron abundance (([Fe/H]) v = 0.01), [Mg/Fe] 
overabundance (([Mg/Fe])v = 0.17) and super-solar 
total metallicity (([Z])v = 0.18 dex) are achieved for 
Salpeter IMF with a star formation episode lasting 
1 Gyr. 

The above example shows that a Salpeter IMF is com- 
patible with the observations provided that tsf is shorter 
than ~ 1 Gyr. Longer tsf lead to overabundances that seem 
too small. 

A flatter slope for the IMF obviously relieves the con- 
straint on tsf. The case of x — 0.7 is shown in the right-hand 
panels of Fig. ||. One can see that tsf as long as 4 Gyr still 
yields an overabundance around 0.29 independent of Z m . 
However, this slope tends also to yield very large [Fe/H] val- 
ues, which do not seem implied by the observations. Even 
with Z m = the average iron abundance is super-solar. 
It should be noticed that this is the consequence of the ex- 
tremely large Zm for this slope, i.e. ~ 17 Zq. However, there 
are no known examples of stars with such high metallicities 
in any environment. 
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Figure 3. Top-left panel: The same as Fig. |2j, however TNH96 nucleosynthesis is adopted. As a consequence, significantly longer star 
formation time-scales 1 Gyr) match the observed range of a-enhancement (0.2 — 0.4 dex). Top-right panel: A projection of the 
contours indicating constant ([Mg/Fe])v to the rgp-x-plane. Bottom panels: Contours of constant iron abundance (left-hand panel) and 
total metallicity (right-hand panel) are plotted as functions of tsf and %■ The contours are not independent of tsf because finite stellar 
lifetimes are taken into account (Schaller et al. 1992). This effect is larger for iron due to the substantial contribution of SN la to the Fe 
enrichment. Generally, the total metallicity is larger than the iron abundance by 0.1 — 0.2 dex. The global values of both Fe/H and Z 
turn out to be sub-solar (—0.4 and —0.2 dex, respectively), if an [Mg/Fe] overabundance is accomplished with Salpeter IMF (x = 1.35). 



Table 5 Moan stpllar ahiindanrng (V-lnminngity weighted) fnr 



the fast clumpy collapse model with tsf = 1 Gyr, x = 0.7 and 
Zm = 0, TNH96 yields. The fraction of gaseous mass transformed 
into stars is varied from /trans = 1.0 (complete exhaustion) to 
/trans = 0.5 to demonstrate the effect on the abundance ratios. 
The fifth column gives the total metallicity of the most metal-rich 
stars formed. 



/trans 


<[Mg/Fc]) v 


<[Fe/H]) v 


([Z])V 


Zm 


1.0 


0.35 


0.45 


0.70 


16.9 


0.9 


0.35 


0.36 


0.62 


12.0 


0.8 


0.36 


0.29 


0.56 


9.9 


0.7 


0.36 


0.19 


0.48 


7.9 


0.6 


0.36 


0.12 


0.41 


6.6 


0.5 


0.36 


0.00 


0.30 


5.0 



This 'overproduction' of metals can be compensated 
by transforming a smaller fraction /trans into stars. Ta- 
ble [| gives the results for various /trans, and fixed r = 1 
Gyr, x = 0.7, Z m = . For this slope, in o r der to obtain 



([Fe/H]>, 



0.2 dex (Worthey et al. 1992 



Davies et al 



199; ), less than approximately 70 per cent of the total 
gaseous mass should be converted into stars, with the rest 
being lost to the intergalactic or intracluster medium. Note 
that the most metal-rich population exhibits metallicities of 
17 Zq for the case of complete gas exhaustion. 



3.2.1 Scalo-like IMF 

Assuming an IMF which flattens out for m 0.5 Mq 
(Scalo 1986; Kroupa, Tout & Gilmore 1993; Gould, Bah- 
call & Flynn 1997) pushes more mass to massive stars and 
thus increases the efficiency of the chemical evolution pro- 
cessing. To investigate the effect on the results we performed 
additional simulations using the estimate of Gould et al. at 
the low-mass end 



m < 0.6 Mq 
0.6 Mq < m < 1 



Mq 



0.44 
-1.21 



and Salpeter slope above 1 Mq. The averaged iron abun- 
dance - as compared to the single slope Salpeter case - is 
then further raised to ([Fe/H])v ~ —0.14 dex, whereas the 
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Figure 4. The figure shows the average stellar [Mg/Fe] overabundance (upper panels) and iron abundance (lower panels) in the fast 
clumpy collapse model as a function of star formation time-scale for Salpeter IMF (right-hand panel) and fiat IMF (left-hand panel). 
TNH96 SN II yields are adopted. The various line-styles denote different minimum metallicities Z m of the central stellar population. 
The symbols represent models excluding Type la supernovae for the two bracketing Z m = (diamonds) and Z m = Zm (triangles). The 
numbers are additionally shown in Table U (columns 4 — 7) . 



Mg/Fe ratio is not affected. Cutting the metal-rich part of 
the CSP at Z m = 1 Zq leads to significantly super-solar 
abundances in iron (([Fe/H])v ~ 0.18 dex) and total metal- 
licity ({[Z]) v * 0.34 dex). 

If instea d a steeper IMF x = 1.7 above 1 Mq is adopted 
(3calo 1986), the iron abundance is slightly increased to 
([Fc/H])v = —0.48 dex. The degree of cv-enhancement how- 
ever, is reduced to ([Mg/Fe])v = 0.1 dex. A Scalo slope of 
x = 1.7 at the high-mass end requires SFTs of the order 10 s 
yr in order to produce ( [Mg/Fe] )v = 0.2 dex. The reason for 
this pattern lies in the fact that apart from the contribution 
of SN la to the Fe enrichment, only high-mass stars above 
8 Mq are responsible for the Mg/Fe ratio. 



3.2.2 The SNIa rate 

The time-scales derived up to now are al so dependent on the 
adop ted description of the SN la rate (Greggio & Renzini 



1985 ). In our models, each stellar generation gives rise to the 
first SN la event after 30 Myr from its birth, and 50% of the 
explosions occur within the first 0.5 Gyr approximately. We 
take this time as a characteristic time-scale for the Fe enrich- 



ment of the ISM from Type la explosions (rFe.ia). The level 
of the SN la rate has been calibrated on the chemical evolu- 
tion of the solar neighborhood (Paper I). In our modeling, 
it yields a current rate of i?i a « 0.05 - 0.1 SNuQ depending 
on the IMF slope and the SFT. This result is in agreement 
with the lowest observational estimates of SNe la rates in 
ellipticals (see Cappellaro 1996, and references therein). 

Given the uncertainties on the SN la progenitors, sub- 
stantially longer TF C ,ia could apply, and correspondingly an 
[Mg/Fe] overabundance could be accommodated with val- 
ues of tsf systematically larger than derived in the previous 
section. We notice however that current models pre dict val- 
ues o f TFc.ia not significantly longer than ~ 1 Gyr (|Greggio 



199C ) . It is instructive to evaluate the possibility of suppress- 
ing completely the Fe contribution from SN la, since it also 
characterizes the assumption of selective mass loss (see op- 
tions ii and iv in the Introduction) . The results are shown in 
Fig. [| as open symbols, where the diamonds refer to the av- 



t lSNu= lSNex(10 10 L o l- 1 X (100 yr) -1 . The luminosities are 



taken from Worthey (1994) 
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eraged abundances of the total metallicity distribution, and 
the triangles to the abundances of the highest Z stellar pop- 
ulation. Note that from stellar abundance ratios it is not 
possible to distinguish the options selective mass loss and 
lower SN la rate. For cluster ellipticals the ICM abundance 
in combination with the stellar abundances could be used to 
distinguish between the two options, under the assumption 
that gEs dominantly contribute to the ICM enrichment. 

Without the Fe enrichment from SN la the {[Mg/Fe])v 
saturates at the maximum SN IFSSP value of 0.28 dex 
(Salpeter) and of 0.4 dex (a; = 0.7), independent of the lower 
Z cut and the SFT, while ([Fe/H])v is shifted to lower val- 
ues by roughly 0.2 dex. Total metallicity (not shown in the 
figure), instead, is lowered merely by 0.01 dex according to 
the small influence of SN la on Z for the SFTs considered 
here. 

Both [Mg/Fe] overabundance and iron abundance be- 
come virtually independent of tsf, since there is no de- 
layed iron enrichment from SN la events. This computa- 
tion trivially shows that without the SN la contribution, 
the a-elements overabundance does not constrain the SFT, 
while it constrains the SSP yields, and therefore, having cho- 
sen the nucleosynthesis prescriptions, the IMF slope. Adopt- 
ing different models for the SN la progenitors with longer 
TFc.ia, and/or assuming that a fraction of the Fe ejected by 
SN la is not incorporated in the ISM (i.e. selective mass 
loss) would lead to intermediate results. Accordingly, longer 
tsf, or steeper IMF slopes would meet the observational 
constraints. A quantitative exploration of the predictions of 
different SN la rates, with the appropriate calibrations on 
the solar neighborhood data and on the observed current 
rates, will be the subject of a forthcoming paper. 



3.3 Central abundances in the model of merging 
spirals 

As a second possibility for elliptical galaxy formation we 
consider the merger of two spirals, and discuss the chemical 
outcome from a star burst triggered by the merging event. 
The main difference with respect to the fast clumpy collapse 
case lies in the initial abundances adopted for the compu- 
tation of the chemical evolution. Our model for the solar 
neighborhood (Paper I) predicts that a half solar and a solar 
metallicity are reached after 3 and 10 Gyr, respectively. We 
take the ISM abundances at these ages as initial conditions 
for the chemical evolution during the burst (see Table |l|). 
We recall here that the model in Paper I adopts Salpeter 
IMF, which allows to reproduce the abundance patterns of 
the solar neighbourhood. 

Modeling the burst of star formation we discuss the 
burst time-scale tsf, however the overall SFT of the ob- 
ject includes both tsf and the age of the parent spirals (see 
Fig. [jj). For consistency, we look at the properties of the final 
galaxy at a total age of 15 Gyr, i.e. when the burst popu- 
lations are 5 Gyr and 12 Gyr old in the cases of solar and 
half-solar pre-enrichment, respectively. Note, however, that 
the passive evolution of the SPs beyond 1 Gyr only provides 
a second order effect on the average abundances. The mini- 
mum metallicities of the populations formed at merging are 
by construction Z m =0.5 Zq and 1 Zq, and represent the 
stars in the galaxy core. The fractional amount of gaseous 



mass is 50 per cent and 20 per cent for the Z[ n = 0.5 Zq 
and Zi n — 1 Zq case, respectively. 

The results of the calculations for the different param- 
eters considered are summarized in the last two columns of 
Table ^, and plotted in Fig. |j| for Salpeter (left-hand panel) 
and flat IMF (right-hand panel). Again, the calculations 
neglecting SNe la are denoted by the symbols (diamonds: 
Z in = 1 Z Q and triangles: Z in = 0.5 Zq). For Salpeter IMF 
both models (solar and half-solar pre-enrichment) fail to pro- 
duce at the same time high metallicity and [Mg/Fe] over- 
abundance. For half-solar initial abundances (characterized 
by an initial a-enhancement of 0.1 dex) the resulting average 
iron abundance is sub-solar. Starting with solar abundances, 
instead, super-solar iron abundance (([Fe/H])v ~ 0.14 dex) 
and metallicity ({[Z])v ~ 0.20 dex) are achieved, whereas 
the [Mg/Fe] overabundance (( [Mg/Fe] )v 5a 0.1 dex) is low. 
The apparently weak influence of neglecting SN la (symbols) 
comes from the fact that in this model SNe la are excluded 
only during the short burst phase. The iron provided by 
SN la explosions during the evolution of the parent spirals 
is included by construction in the initial conditions. Clearly, 
a flatter IMF during the burst yields a large overabundance, 
and high metallicities. It is worth noticing that the duration 
of the burst is not important with respect to the constraints 
considered here. 

3.4 Comparison between the two models 

We now turn to compare the results in the two extreme 
modes for the formation of elliptical galaxies. Inspection of 
Table reveals that, for the same degree of pre-enrichment 
(i.e. Z m ): 

(i) The average metallicity ([Z])v of the central CSPs is 
approximately the same for the two models. 

(ii) The ([Fe/H])v is systematically larger, and the aver- 
age [Mg/Fe] overabundance is systematically lower, in the 
merging spirals model. 

The first result simply reflects the fact that in both cases 
we have an overall closed box metallicity distribution with 
the same widths, the maximum metallicity being determined 
only by the nucleosynthesis prescriptions and the IMF slope. 
The second result is a consequence of the different initial 
conditions: in the merging spiral model all the Fe ejected 
by SN la during the past several Gyr of evolution of the 
parent spirals is incorporated in the gas undergoing the fi- 
nal burst of SF in the centre. In the fast clumpy collapse 
model, even the high metallicity stars form early, when pol- 
lution by SN la is still relatively modest. The larger a- 
enhancement achieved in this model is a direct consequence 
of the same effect. Indeed, the differences in both ([Fe/H])v 
and ( [Mg/Fe] )v are larger for the CSPs with Z m = 1 Zq, 
with respect to the cases with Z m = 0.5 Zq. When tsf is 
sufficiently long (e.g. 4 Gyr) the two models converge to sim- 
ilar results for ([Fe/H])v and ([Mg/Fe])v, since in this case 
also in the fast clumpy collapse model most of the Fe ejected 
by SN la has been included in the high metallicity tail of the 
CSP. 



4 DISCUSSION 
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Figure 5. The figure shows the average stellar [Mg/Fe] overabundance (upper panels) and iron abundance (lower panels) in the merging 
spirals model as a function of burst time-scale for Salpeter IMF (right-hand column) and fiat IMF (left-hand column). TNH96 SN II 
yields are adopted. The line-styles are the same as Fig. ^, denoting the respective degrees of pre-enrichment (see also Table [l]). The 
symbols represent models excluding Type la supernovae during the burst for Z m = 0.5 Zq (triangles) and Z m = 1 Zq (diamonds). The 
numbers are additionally shown in Table U (columns 8 — 10). 



4.1 Fast clumpy collapse vs. merging spirals 

The two formation scenarios considered here can be distin- 
guished by the time-scale on which the pre-enrichment of the 
central SPs to solar-like abundances is accomplished. In the 
merging spirals scenario, this time-scale is as long as several 
Gyr, and the pre-enriched gas that forms new stars has so- 
lar element abundances, in particular non a-enhanced abun- 
dance ratios. In the fast clumpy collapse scenario, instead, 
the pre-enriched gas forming the nucleus is overabundant 
in a-elements because of the preceding short SFT (shorter 
than 1 Gyr). 

Our computations show that for the clumpy collapse 
both constraints for gEs, namely a-enhancement and super- 
solar metallicity, can be matched under the assumption of 
Salpeter IMF and SFTs up to 1 Gyr. High average metal- 
licities are accomplished through segregating the metal-rich 
populations (with a minimum metallicity around 1 Zq) in 
the central region from the rest of the galaxy. 

Super-solar iron abundance can be obtained if in ad- 
dition a flattening of the IMF at the low-mass end around 
0.5 Mq and Salpeter slope for masses above 1 Mq are ap- 



plied. Flatter IMFs (in particular at the high mass end) 
generally allow longer SFTs, for the same ([Mg/Fe] )v, and 
lead to larger average metallicities. The metal 'overproduc- 
tion' can be avoided by assuming that only a fraction of the 
initial gas mass is converted into stars, with consequences 
on the chemical enrichment of the intergalactic (or intra- 
cluster) medium. Besides, flatter IMFs, with metallicity dis- 
tributions extending up to several times Zq, do not require 
the segregation of the high metallicity SPs in the galaxy 
cores, to accomplish average metallicities in agreement with 
the observational indications. It is worth noticing, however, 
that there is no guarantee that these metallicity distribu- 
tions would yield the strong Mg2 indices observed in the 
central parts of the galaxies. In fact, the value of the indices 
depend on the fitting functions which seem to saturate at 
large metallicity (see G97). 

In the case of merging of two Milky Way galaxies, the 
only way to achieve a high a-enhancement is to invoke a flat 
IMF in the final burst. This implies a large average metal- 
licity for the central SP, with {[ZT])v ~ 0.5 already for x=l. 
In this case, it seems unlikely that the final burst does not 
proceed up to gas exhaustion, since the chemical process- 



© 0000 RAS, MNRAS 000, 000-000 



a- enhancement in elliptical galaxies 11 



ing takes place deep in the potential well. Thus, if ellipticals 
form in this way, IMF slopes flatter than x = 1 seem highly 
unfavoured from the values of the Mg and Fe indices ob- 
served. The only way out would be to assume very short 
time-scales for the final burst of Tburst ~ 10 7 yr, which ap- 
pears contrived. The merging of two spirals at early times 
in their chemical evolution looks more favourable, since our 
models starting from half solar initial conditions result in 
abundances marginally consistent with both observational 
constraints (a-enhancement and high metallicity) . 



4.2 Gradients inside the galaxy 

In the fast clumpy collapse model, in order to accomplish 
large average metallicities in the core population we assumed 
a cut in the Z distribution, such that only stars with Z > Z m 
are found in the central parts of the galaxy. Correspondingly 
the low metallicity stars (Fig. ^J) are in the outskirts, which 
qualitatively accounts for the observed metallicity gradients. 
As a consequence one expects that the stars in the outskirts 
of the elliptical galaxy should be more a-enhanced, with 
respect to those in the central parts, i.e. we expect a positive 
gradient in the [Mg/Fe] ratio. This pattern is present in the 
solar neighbourhood, with the metal-poor stars being more 
a-enhanced. In the merging spirals model, instead, the stars 
in the outer parts of the galaxy have approximately solar 
abundance ratios, i.e. ellipticals formed in this mode should 
exhibit negative gradients of [Mg/Fe]. 

Although the derivation of real abundances from line 
strengths measurements is highly uncertain (G97; Tantalo, 
Bressan & Chiosi 1997a, 1997b), so far the observations tend 
to indicate that the a-enhancement is either constant, or 
mildly dec reasing from the centre to the ou ter parts of the 
ellipticals ( pavies et al. 199$ |Mehlert 199g| ). A trend of in- 
creasing overabundance with increasing metallicity (with a 
constant IMF slope) can only be achieved in a scenario in 
which the formation of the metal-rich population is evolu- 
tionary decoupled from the star formation episode which 
generated the low metallicity stars, as assumed in our model 
of the merging spirals. In principle, one could apply this 
same decoupling also to the fast clumpy collapse picture, 
assuming that the SN la products of the low Z stellar gener- 
ation do not mix efficiently with the collapsing gas, and are 
preferentially lost to the IGM. However, this looks rather 
artificial, due to the very short time-scale over which the 
formation of the whole galaxy is accomplished. Anyway, the 
measurement of abundance gradients offers the opportunity 
to discriminate between these two formation scenarios. It is 
therefore very important to better understand how the line 
strengths map into real abundances, and to study the trend 
of [Mg/Fe] overabundance with radius in great detail and in 
many objects. 



are adopted from Worthey ( |l994| ). In agreement with esti- 
mates by Edmunds (1992), the metallicity averaged by mass 



turns out to be systematically higher by 0.1 dex for Salpeter 
IMF, independent of the SFT: 

A[Fe/H] = ([Fe/H])» - ([Fe/H])v ~ 0.1 dex 



4.3 Caveats and limitations 

4-3. 1 Luminosity vs. mass average 

It is usually argued that the diffe rence between luminos- 
ity and mass av erage is negligible ( Arimoto fe Yoshii 1987 
Matteucci 1994 ). The derived abundances given in this pa- 
per are all weighted by the light emitted in the V-band, the 
specific V-luminosities as functions of age and metallicity 



Following the argument by Arimoto & Yoshii (1987) this is 
due to metal-poor K giants dominating the light in the V- 
band (see also Edmunds 1990). The effect on the [Mg/Fe] 
overabundance is significantly smaller, i.e.: 

A [Mg/Fe] = ([Mg/Fe]), - ([Mg/Fe]) v ~ -0.01 dex 

A flattening of the IMF lessens the discrepancy between 
mass and ligh t ave rage, hence A[Fe/H] ~ 0.05 dex. How- 
ever Worthey (1994) provides models for Salpeter slope only, 
causing an inconsistency in our approach. Indeed, at large 
ages the stellar M/Lb ra tios are affected by the IMF slope 
adopted (Maraston 1995). To check the influence on our re- 
sults, we additionally perfor med s imulations using SSP lumi- 
nosities derived by Buzzoni (198£ ), who provides models also 
for flatter IMF slope. It turns out that a A[Fe/H] « 0.05 dex 
for Salpeter IMF is further reduced to A[Fe/H] at 0.02 dex 
with x = 0.7. We conclude that within an uncertainty of 0.1 
dex for the metallicity the difference between mass and light 
average in fact provides a minor effect independent of the 
IMF slope. 



4.3.2 Stellar Yields 

Our conclusions are strongly affected by uncertainties of 
SN II nucleosynthesis (Paper I; Gibson 1997a, 1997b; Porti- 
nari, Chiosi & Bressan 1998). We analyse the consequence 
of this problem by initially considering two different SN II 
yields (WW95, TNH96), but most of the computations are 
performed with TNH96, since they give the better fit to the 
solar neighbourhood data (Paper I). However, since TNH96 
underestimate the Mg/O ratio by roughly 10 per cent (Pa- 
per I), and since the oxygen yield is likely to be better de- 
termined, there is the possibility that the Mg stellar yield 
is underestimated in TNH96 models. This leaves room to 
further relax the constraints on tsf and /or to reconcile the 
a-enhancement in ellipticals with the local IMF. 

In our computations we assume an upper mass cut-off 
of 40 M Q for both WW95 and TNH96 SN II yields. In- 
creasing the upper mass cut-off to 70 Mq , and interpolating 
betwee n TN H96 yields for a 40 Mq model and the Nomoto 
et al. (1997) yields of a 70 M0-star, leads to an increase 
of 



0.15 dex of the [Mg/Fe]-ratio in the SN II ejecta in 
the case of Salpeter IMF (Paper I). Again, this relieves the 
constraints on the short tsf for the clumpy collapse model, 
or on the IMF slope. In particular, for the case of solar 
pre-enriched gas (merging Milky Ways) we performed addi- 
tional calculations including the 70 Mq contribution to the 
ISM enrichment. Indeed, for Salpeter IMF we now obtain 
( [Mg/Fe] )v ~ 0.2 dex for tsf = 0.3 Gyr . However, the ac- 
tual contribution to the ISM from stars more massive than 
40 Mq is highly uncertain (Paper I), due to the fall-back 
effect (WW95). Anyway, if these stars contribute substan- 
tially to the ISM pollution, the discrepancy in the derived 
overabundance between the fast collapse and merging, spirals 
models becomes larger than what is shown in Table Thus, 



© 0000 RAS, MNRAS 000, 000-000 



12 D. Thomas, L. Greggio and R. Bender 



the conclusion that these two extreme cases of galaxy for- 
mation scenarios lead to different abundance ratios is even 
reinforced. 



5 CONCLUSION 

We have explored how to achieve the observed a- 
enhancement in luminous ellipticals in the three-dimensional 
parameter space of star formation time-scale, IMF slope, and 
stellar nucleosynthesis. As anticipated in Paper I on the ba- 
sis of SSP yields, we find that, using the stellar yields from 
WW95, the observed {[Mg/Fe])v are reproduced only with 
extremely short star formation time-scales (SFT) of the or- 
der 10 7 yr. Therefore most of the results presented here have 
been obtained using TNH96 models. 

It turns out that the abundance ratio ([Fe/H])v under- 
estimates the real total metallicity of the system by 0.1 — 0.3 
dex depending on the SFT and IMF slope. We also briefly 
explored the difference between (V-band) light and mass 
average metallicity, finding that the mass averaged [Fe/H] 
abundance ratios are w 0.1 dex larger than the V-luminosity 
averaged ones. The effect on [Mg/Fe] is instead negligible. 

We have considered to extreme modes of elliptical 
galaxy formation, the fast clumpy collapse and merging spi- 
rals modes. We find that in the fast clumpy collapse model 
( [Mg/Fe] ) v ~ 0.2 dex (in the high metallicity SPs) is ob- 
tained with IMF slopes close to Salpeter, provided that the 
overall formation time-scale does not exceed 1 Gyr. For the 
merging spirals model, instead, a sizeable overabundance is 
achieved only by claiming a flat IMF for the burst forming 
the central population. This is due to the large Fe content of 
the ISM at the epoch of the burst. The earlier the merging 
takes place in the chemical evolutionary history of the par- 
ent spirals, the easier it is to accomplish the overabundance 
in the burst population. In general, galaxy formation time- 
scales exceeding a few Gyr do not provide ([a/Fe])v ratios 
consistent with observations, unless a significant deviation 
from the local IMF is allowed. 

The two extreme models tend to predict opposite trends 
for the a-enhancement gradient within the galaxy: the 
fast clumpy collapse leads preferentially to larger overabun- 
dances in the outskirts where the low Z stars are formed. 
For the merging spirals model we expect approximately solar 
ratios in the outer stellar populations. If the merging occurs 
relatively early in the chemical evolution of the parent spi- 
rals, a slight a-enhancement is expected also in the outer 
stellar populations. Thus the behaviour of { [Mg/Fe] )v with 
radius can give imp ortant hints on the galaxy formation 
process. Kauffmann ( 199€j ) proposed an intrinsic difference 
between the formation of cluster and field ellipticals, due 
to the high density environment favouring shorter galaxy 
formation (total) time-scales. If field objects are formed in a 
later merging event with respect to cluster objects we expect 
that field ellipticals tend to exhibit lower a-enhancement 
and steeper a-enhancement gradients than cluster ellipti- 
cals. 
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